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Abstract 

We prove the existence of discrete breathers (time-periodic, spatially localized solu- 
tions) in weakly coupled ferromagnetic spin chains with easy-axis anisotropy. Using 
numerical methods we then investigate the continuation of discrete breather solutions 
as the intersite coupling is increased. We find a band of frequencies for which the 1- 
site breather continues all the way to the soliton solution in the continuum. There is 
a second band, which abuts the first, in which the 1-site breather does not continue 
to the soliton solution, but a certain multi-site breather does. This banded structure 
continues, so that in each band there is a particular multi-site breather which continues 
to the soliton solution. A detailed analysis is presented, including an exposition of how 
the bifurcation pattern changes as a band is crossed. The linear stability of breathers 
is analyzed. It is proved that 1-site breathers are stable at small coupling, provided a 
non-resonance condition holds, and an extensive numerical stability analysis of 1-site 
and multisite breathers is performed. The results show alternating bands of stability 
and instability as the coupling increases. 
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1 Introduction 



Discrete breathers are time periodic, spatially localized solutions in networks of coupled 
oscillators (including rotors and spins). They arise in a variety of very general systems 
13, 16 1 due to the interplay between nonlinear and discrete effects, and therefore there is 



an enormous potential for their application in many areas of physics, particularly condensed 
matter and biophysics, where physical systems are often inherently discrete. Generally one 
finds that as the system is moved closer to its continuum limit, by increasing a coupling 
constant in the theory, there comes a point at which a discrete breather solution no longer 
exists. This is to be expected since an increase in the coupling constant results in an 
expansion of the phonon frequency band and eventually this band captures a harmonic of 
the breather frequency. In fact, it is often the case that the discrete breather cannot be 
continued even up to the point at which the above resonance argument applies, and this is 
not well understood. 

In this paper we study the classical equations of motion for a ferromagnetic spin chain 
with an easy-axis anisotropy. Some numerical studies of discrete breathers (which are known 
as intrinsic localized modes in the condensed matter physics literature) have been performed 
in the cases of easy-plane ferromagnets |17] and easy axis antiferromagnets |12]. However, 
both the perspective and the results of the current investigation are quite different. We 
prove the existence of discrete breathers in the weakly coupled case by starting from the 
anti-continuum limit (zero coupling constant) and applying an implicit function theorem 
argument to the trivial 1-site breather. We then apply numerical methods to investigate 
the continuation of this solution as the coupling constant is increased, with the novel result 
that it continues for all values of the coupling. This is explained by an examination of 
the spin waves. Depending upon the breather frequency, the continuation results in either 
the soliton solution of the continuum model or a trivial static solution in which all spins 
are either aligned or anti-aligned with the vacuum. Similar results are found to apply to 
multi-site breathers, which fit into an intricate pattern, and lead us to conjecture that 
the continuum soliton solution (of any allowed frequency) can always be obtained from 
the continuation of a particular form of multi-site breather. Evidence in support of this 
conjecture is presented. The numerical results are in some ways reminiscent of those found 
[^, ^] for the discrete nonlinear Schrodinger (DNLS) system, but have not previously been 
seen in spin chains. 

We then go on to perform a linear stability analysis of the breathers within the Krein 
theory framework pioneered by Aubry iQ. We prove that nonresonant continued 1-site 
breathers must be stable for sufficiently weak coupling, and numerically investigate the 
stability of both 1-site and multisite breathers. The results suggest that those families 
which tend to the continuum soliton experience regular repeating bands of instability of 
diminishing strength as the coupling increases. 

During the refereeing process for this paper, a very interesting paper by Flach, Zolotaryuk 
and Fleurov on breathers in very similar systems appeared Q. The differences between 
their results and ours are, briefly, as follows. They consider both easy-axis and easy-plane 
anisotropy, but in the case where the exchange interaction is anisotropic also (our exchange 
interaction maintains isotropy). Breather existence is proved by continuing from a limit in 
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which two of the exchange integrals vanish, but the third remains nonzero. Since the contin- 
uation is for small values of the continuation parameters (the two small exchange integrals) 
their existence result applies to strongly anisotropic exchange interaction, as opposed to 
ours, which applies to isotropic exchange. More importantly, in this paper we exploit the 
isotropic exchange to reduce the breather equations to a purely algebraic system, which 
greatly simplifies the analysis in comparison with theirs. Consequently our existence result 
requires none of the nonresonance hypotheses of Flach et al's. Also we are able to build spa- 
tial localization of breathers directly into our analysis, whereas Flach et al do not consider 
this issue at all (although standard results of MacKay on exponential breather localization 
should apply to their system just as they do to ours - see section p. Similarly, we exploit 
the exchange isotropy to simplify the linear stability analysis, again reducing it to a purely 
algebraic problem, a reduction not possible for the systems analyzed in |^]. We are able, 
therefore, to perform a very detailed numerical existence and stability analysis at minimal 
computational cost. To summarize, Flach et al prove a result of rather more general appli- 
cability and physical relevance than ours, though their context does not, strictly speaking, 
include ours. We consider a somewhat more idealized system, but obtain correspondingly 
stronger results. 

2 Anisotropic ferromagnetic spin chains 

The classical formulation of a spin chain involves a three-component unit vector rij, 
giving the spin at each lattice site i G Z. The type of spin chain is defined by its Hamiltonian, 
which we take to be 



Here a > is the coupling constant (exchange integral), which is positive in the case of a 
ferromagnet, and ^ > is the anisotropy constant, which is also positive since we wish to 
consider an easy-axis anisotropy, with = (0, 0, 1) being the easy-axis. In this case the 
minimum of the Hamiltonian is zero, obtained by the vacuum configuration m = ±63. In 
the following we shall choose -1-63 to be the vacuum configuration and refer to spins which 
take the values -1-63 and —e^ as being spin up and down respectively. Furthermore, by 
applying a scaling symmetry we can, without loss of generality, set A = 1. 
The equation of motion is obtained from the Hamiltonian as 




(2.1) 



OH 

Ui = -Ui X -— 



(2.2) 



where a dot denotes differentiation with respect to time. Using ( |2.1| ) we obtain 




(2.3) 



Discrete breather solutions have the form 




(2.4) 
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where uj is the frequency. Substituting this ansatz into the equation of motion ( |2.3D yields 
the following nonlinear second order difference equation for the angles 9i, 

a{cos9i(sm6i^i + sin6'j_i) — sin 0.j(cos ^i+i + cos0j_i)} = sin0j(cos6'j — u). (2.5) 

Such a periodic solution may properly be called a discrete breather if it is spatially localized, 
that is, limj^-i-oo = 0. 

In section H we shall prove the existence of discrete breathers for a sufficiently small 
and in section ^ we shall study them numerically for a range of values of a. 

Note that for all a there is a set of static solutions in which each spin can independently 
be chosen to point either up or down, since in each case we have that sin^j = for all i, 
which clearly solves equation ( p. 51 ). These solutions will play an important role later in our 
discussion. 

Finally, in this section, we address the continuum limit a oo. Write a = l/h? and 
regard iij as the value of a continuous function n(a;) sampled at the lattice points x = ih. 
Taking the continuum limit /i ^ 0, the Hamiltonian (2.1) becomes 



H = |{l„'.n' + i[l-(n,.e3)2]}dx (2.6) 
and the equation of motion is 

h = n X n" + (n • e3)(n X ea). (2.7) 

where prime denotes differentiation with respect to x. 

This partial differential equation has time periodic, exponentially localized solutions 



known as magnetic solitons [11|. At this point it is perhaps worth pointing out that there is 
an unfortunate difference in nomenclature, in that in the discrete case time periodic localized 
solutions are known as breathers, whereas in the continuum model they are termed solitons. 
We shall continue to use the term soliton when referring to the continuum limit, but the 
reader should be aware that it has exactly the same time dependence as the discrete breather 
and moreover, as we shall see later, the soliton can be obtained from the discrete breather 
in the continuum limit. 



Using the same form for the time dependence as in the discrete case (|2.4| ) i.e. 

n{x,t) = (sin 0(x) coswi, — sin6'(x) sinwt, cos 6'(x)) (2-8) 
equation ( p. 7] ) yields the second order ordinary differential equation 

6»" = (cos 0-w) sine (2.9) 



which, of course, is also obtained from the continuum limit of ( |2.5| ). The boundary condi- 
tions for a soliton solution, located at the origin, are ^'(0) = 6{oo) = 0. Equation ( p.S| ) can 
be integrated explicitly and the solution satisfying the correct boundary conditions is given 
by 

e{x) = cos-H— .^'',2, . - 1}. (2.10) 

1 — (1 — uj) tanh (xvl — "^j 
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Since we require the soliton to be a smooth exponentially localized solution, this formula 
shows that the frequency must be restricted to the range a; E (0, 1). There is a second type 
of soliton solution [0] for which a; < 0, but this has a different structure from ( |2.10 ); in 
particular n = —e^ at the centre of the soliton for all uj. This negative frequency soliton 
will not arise in our discussion, so in this paper when we refer to a soliton we shall mean 
the solution ( |2T0| ). 



3 Analytic results on discrete breathers 

Following the "homoclinic orbit" approach of Flach one would like to prove directly 
the existence of solutions of ( |2.5| ) with the correct boundary behaviour, using techniques of 
dynamical systems theory. Unfortunately, equation ( |2.5D does not determine a well-defined 
homeomorphism of the torus, {6i-i, 9i) i— > {6i, 9i+i), so the direct approach is not convenient 
here. 

However, we may still prove the existence of breathers in this system by continuation of 1- 
site breathers from the decoupled limit (a = 0), in the manner of Aubry and MacKay's work 
on oscillator networks [ll3|| . Here we have an algebraic system rather than an infinite system 
of ODEs so the details are considerably less technical (compare, for example, with |^]), and 
will be treated with corresponding brevity. The idea is that when a = 0, a; G (— 1, 1), ( ^.5| ) 
supports the almost trivial solution 



i / 

cos^"*^ Lu i = 



-1,, ._n (3-1) 



in which every spin remains pointing up except one (whose location we have chosen to be 
1 = 0) which precesses with frequency uj around a circle of fixed latitude. Note that there 
is no reason to assume that w > as we had to in the continuum system. Keeping u fixed, 
the existence of breathers for all a sufficiently small will follow from an implicit function 
theorem argument. 

To be precise, let F : £2 ®M ^ £2 such that 

[F{9,a)]i = (cos — w) sin^j — a[cos 6'j(sin ^j+i + sin6'j_i) — sin6'j(cos ^.j+i + COS (3.2) 

where £2 is the Banach space of sequences 9 : Z ^M. with finite norm 



(3.3) 



Note that ^ G ^2 is a solution of ( ^.5]) at coupling a if and only if F{9, a) = 0. Note also 
that every 9 £ £2 converges to as i — > ±00, so such a zero of F(-, a) is a discrete breather. 
For all to G (—1,1), F{9^,0) = 0, where 9^ is the 1-site breather defined in (3.1), and the 
partial differential of F with respect to the £2 factor at {9, a) = {9^, 0) is easily seen to be 
an isomorphism DFf^Qb Q-^ : £2 £2, that is, a bijection with bounded inverse. In fact this 
linear map is diagonal: 

Hence we may apply the implicit function theorem P] to obtain the following 
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Theorem 1 For all uj G (—1,1) there exist e > and a map [0, e) 
a 1-^ 0", such that 9°^ is a frequency to, coupling a solution of (2^) and 9^ 
9°^ is unique provided e is chosen sufficiently small. 



£2, denoted 
9^. The map 



As noted above, since the continuation occurs in £2, it has weak spatial locahzation built 
in. This can be improved to exponential spatial localization (that is, there exist C, A > 1 
such that \9i\ < CA"'*') by applying some standard results of Baesens and MacKay. The 
idea is to differentiate the continuation equation F{9°^,a) = with respect to a to obtain 
an ODE on £2, namely, 

— = [I)F(,.,„)]-iA(r), (3.5) 



where /\ : £2 ^ £2 such that 

[A(6')]i = cos 6'j(sin 6'j_|_i + sin 6'j_i) — sin ( 



(cos i 



1 + cos ( 



(3.6) 



The theorem above may be interpreted as asserting the local existence and uniqueness of a 
solution to the initial value problem 9^ = 9^ for ( |3.5| ). Now /S.{9°^) is clearly exponentially 
localized provided 9" is, as is q)]^^A(0'^) by theorem 3 of |^]. Thus one obtains an 

exponential decay estimate for d9°'/da at ao given a decay estimate for 9°"^. This may be 
integrated to show that the estimate for 9°"^, if chosen correctly, remains valid in a finite 
interval containing ao (see the proof of theorem 2 of Q , for example) . It remains then to 
note that the initial datum 9'' trivially satisfies all exponential decay criteria (i.e. one can 
choose any A > 1), so exponential localization persists for a sufficiently small. 

We may deduce the optimal exponent A for given lo and a by linearizing equation (^ 
about 9 = 0: 

2 + 



59,. 



a 



69i + 69i 



0. 



(3.7) 



The general solution to (|3.7D is 



A = 1 + 



AX" + BX-^ where 



UJ 



2a 



+ 



LO 



+ 



a 



■ UJ 



2a 



(3.8) 



A 1*1, which we will encounter in our numerical 



It is this type of exponential decay, \9i 
analysis. 

Another interesting piece of information can be deduced from equation (|3.5|), namely 
the tangent vector to the continuation curve in I2, that is, the direction of continuation 
away from a = 0. Substituting (3^) and ( |3.4D into (3^) at a = yields 







da 


a=0. 





l+g; 



\i\ > 1 
\i\ = 1 

i = 0. 



(3.9) 



So the continuation initially proceeds by pulling the central spin, no, further away from the 
vacuum es, and pulling its nearest neighbours, n-ti away from es in the same direction (as 
no and each other), while leaving all other spins fixed at e^, to first order. 
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We have proved the existence of a continuation of the one site breather 0^, but the same 
argument can be apphed to any zero of 0) : ^2 — > ^2 to prove existence of more general 
discrete breathers at small a. Clearly, if F{6,0) = then 6i S vrZ U cos~^{w} for all i, 
and given that 6 E. £2, &i = for all \i\ sufficiently large. So one can construct frequency lv 
solutions at a = like 

6* = (..., 0, 0, cos~^ uj, vr, tt, vr, cos^"*^ uj, 0, 0, . . .) (3.10) 

for example. It is easy to check that DF^^^ is an isomorphism, so an analogous result 
to the theorem above applies here also. Of particular relevance to the numerical work in 
section ^ will be those periodic solutions obtained by continuing a = solutions with even 
reflection symmetry about i = 0, that is, 6-i = 9i. Such a solution may be specified by a 
finite coding sequence |14| representing the values of 6i for i = 0, 1, 2, . . . by + if = 0, 



o \i 9i = cos~^ti; and — if 0, = vr, with the convention that the last symbol represents 
the constant tail of the sequence. Of course, since £ £2, each coding sequence must end 
in +, so the last symbol is redundant. Nevertheless, we include it in order to make the 
notation more suggestive. The symbols themselves are meant to represent the directions of 
the central spins: up is +, down — , while o is somewhere (a; dependent) in between. We 
shall refer to those breathers obtained by continuing such a solution by the same coding 
sequence. Hence 0°^, the continuation of 6^ is a (o+)-breather, while 9 and its continuation 
are ( o+)-breathers, for example. 



4 Numerical results on discrete breathers 

In this section we present the results of an extensive numerical investigation of the 
continuation of 1-site and various multi-site breathers from the anti-continuum limit. To 
construct numerical solutions we must first truncate to a finite number of lattice sites, which 
we implement by fixing vacuum boundary conditions, 6i = for i = iN, and restricting 
to the interior of the lattice |i| < A^. Furthermore, the solutions we consider are symmetric 
about the central lattice site, that is, 6i = 6-i, so only the sites i > need be considered, 
with an appropriate modification at the central site i = 0. Explicitly, the task of finding a 
numerical solution reduces to finding a zero of the following (A^ -|- l)-component vector 

{a[cos9i{sm6i^i + sin6'j_i) — sin0j(cos0i+i -|- cos0j_i)] -|- sm6i{uj — cos 9i) < i < N 
2a[cos ^0 sin 61 — sin Oq cos 61] + sin 9o{uj — cos 60) i = 

On i = N 

(4.1) 

as a function of the {N -|- 1) values 6i, i = 0, .., N. 

To find a zero of the vector Fi, for a given value of a, we begin with the value a = 0, 
where we have the explicit solution corresponding to the trivial 1-site breather ( |3.1D (and 
later we use other multi-site breathers). We then increase a by a small amount and use a 
Newton-Raphson scheme to converge to the new solution. This process is repeated until 
the desired value of a has been obtained. During the calculation the determinant of the 
Jacobian matrix Jij = dFi/dOj is monitored to ensure that it is non-zero, as a check that 
the scheme is not accidentally jumping to a different solution branch. The results presented 
below were obtained from the value N = 20, though different size lattices were checked to 
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confirm that, providing N is sufficiently large, the results are not sensitive to the number 
of sites. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

alpha 



Figure 1: Plots of cos^o and cos^i for the continuation of the (o+)-breather with uj = 0.2, 
from a = to a = 0.5. 

In Figure 1 we display the results of the continuation of the (o+)-breather, with fre- 
quency u! = 0.2, from a = to a = 0.5. Presented are plots of cos^o and cos 6*1, that is, 
the third component of the spin for the central and next to central sites. From this plot 
it is clear that the (o+)-breather continues until it joins the static (— +) solution. This is 
true for all frequencies in the band u; G (—1,0.424], where we have computed the edge of 
the band to within the numerical accuracy given. It is not surprising that for uj close to — 1 
the (o+)-breather continues to the static ( — h) solution since for a = the (o+)-breather 
tends to the ( — h) solution as a; — > — 1. 




-0.4 I ' ' ' ' ' ' ' ' ' 1 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

alpha 

Figure 2: Plots of cos^o and cos^i for the continuation of the (o+)-breather (solid curves) 
with ui = 0.5, from a = to a = 2. The dashed curves are the corresponding quantities for 
the soliton solution. 

In Figure 2 we present the results of a similar calculation, but this time the initial (o+)- 
breather has a frequency u; = 0.5. The values of cos and cos 9i are shown as solid curves 
and the dashed curves represent the same quantities for the continuum soliton solution 
( |2.1C| ) sampled at lattice sites x = and x = h, where h = As the coupling a is 

increased, the values at these sites tend towards those of the continuum soliton solution, 
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Figure 3: cos 9i for all the lattice sites of the lo = 0.5 (o+)-breather with a = 2 (diamonds). 
The dashed curve is the corresponding soliton solution. 



demonstrating that the (o+)-breather approaches the soliton solution in the continuum 
limit. As further evidence, in Figure 3 we plot cos^j for all lattice sites, for the value a = 2, 
superimposed with the soliton solution given in ( ^.lOD , with x = i/\/2. As a final check, 
starting with initial conditions obtained from sampling the soliton solution at the lattice 
sites and continuing backwards to a = 0, we find that we recover the (o+)-breather. We 
have also verified that the decay rate of the (o+)-breather shown in Figure 3 fits extremely 
well to the derived exponential decay constant given in 

Similar results apply for all frequencies in the band lo £ (0.424, 1), with the (o+)- 
breather continuing to the soliton solution, although the region in which uj is close to 1 is 
numerically inaccessible within the present scheme since the soliton has a relatively slow 
spatial decay in this region. 

Just below the interface of the two bands, that is for < 0.424 — a; ^ 1, there is 
a more complicated bifurcation structure between the (o+)-breather and the static ( — h) 
solution than in the bulk of the lower band (—1,0.424]. To examine this region requires a 
slightly more sophisticated numerical approach, as we now describe, by illustrating the case 
oj = 0.42. 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 4: cos^o ^or the continuation of the (o+)-breather (dashed curve), the static ( — h)- 
solution (solid line), and the perturbed ( — |-)-solution (dotted curve). 
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00 = 0.424 



(0 = 0.133 
CO = 0.045 
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Figure 5: A table displaying the results, as uj is varied, of continuing breathers of type 

(o+),(— o+) and ( o+). The symbol at the foot of each column indicates which breather 

is being continued and the symbol inside the column represents the end result of the con- 
tinuation at that frequency. 



CO = 0.600 



CO = 0.425 



CO = 0.420 



CO = 0.200 







Figure 6: A schematic to indicate how the bifurcation pattern between the ( — h) solution 
and the (o+) and (— o+) breathers varies with oj. 



For UJ = 0.42 the continuation of the (o+)-breather is presented in Figure 4, where 
the dashed line represents the value of cos^o- The continuation fails at a = a* = 0.4085, 
where the determinant of the Jacobian matrix J is zero. We expect that this solution 
still bifurcates from the static ( — h) solution, but that it first turns around in a space. 
To confirm this expectation we need to find the value a = a, at which the bifurcation 
of the ( — h) solution takes place, and compute the tangent vector in the direction of this 
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bifurcation. The Jacobian of the ( — h) solution is the tridiagonal matrix 
/2a-uj-l -2a 



J 



—a 



uj-l 
a 



a 



-2a + u- 1 



a 



a 



V 



-2a + u- 1 





a 
1/ 



(4.2) 



Computing the determinant of this matrix (for N = 20) and setting uj = 0.42, we find 
that det J has only one real and positive root, which occurs at a = a = 0.3858. Note that 
this calculation provides the first check on our bifurcation assumption since we require that 
a < a* , which is indeed true. Substituting a = a into J we then compute the eigenvector, 
56, corresponding to the zero eigenvalue, normalized so that its norm \59\'^ = e <C 1. Finally, 
we take the ( — h) solution, in which Oq = tt and 9i = for i > 0, and add the perturbation 
66, to create an initial condition which we then iterate to converge to a solution for a value 
of a obtained by incrementing S by a small amount. Following this solution for increasing 
a wc find that it continues until a = a* , where it meets the end point of the continued 
(o+)-breather. This is illustrated in Figure 4, where we plot cos (dotted curve) for the 
solution branch obtained by perturbation of the static ( — h) solution, as just described. 
Following the end point of this solution branch backwards in a confirms that it joins the 
static (— +) solution at a = q. 

To summarize, the results of the continuation of the (o+)-breather, as lo is varied, are 
represented in the first column of Figure 5, and presented schematically in the upper half 
of Figure 6. For u G (0.424, 1) the (o+)-breather continues to the soliton, for u outside this 
range, but close to 0.424, it turns back and joins the static (— +) solution, and in the bulk 
of the band (—1,0.424] it joins the static (— +) solution without turning back. 



Figure 7: cosi 
a G [0,0.8]. 



, cos ^1, cos ^2 for the continuation of the u> = —0.2 (— o+)-breather for 



For u} above 0.424 there is no longer a bifurcation of the static ( — h) solution with 
the (o+)-breather, but to fill in the full bifurcation pattern of the ( — h) solution we must 
now consider the continuation of the (— o+)-breather. Note that now the start point of 
the continuation, a = 0, is already a multi-site breather, since we have applied reflection 
symmetric conditions, 9i = 9-i, and the breathing site is not located at i = 0. 
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Figure 8: cos cos ^i, cos 6*2 (solid curves) for the continuation of the uo = 0.2 (— o+)- 
breather for a G [0,2]. The dashed curves are the corresponding quantities for the sohton 
solution. 



0.5 - 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 9: cos 9o, cos 9i for the continuation of the u = 0.55 (— o+)-breather for a G [0, 0.5]. 



In Figure 7 we plot cos 9i for the first three sites, i = 0, 1, 2, for the continuation of the 
(—o+) -breather with frequency u = —0.2. Clearly, for this frequency, the (—o+) -breather 

continues to the static ( h) solution, and this is true for all uj G (—1,0.133). From 

a similar earlier discussion, it follows that this result is to be expected, at least for uj 
near —1. For cu G [0.133,0.424] the (—o-f) -breather approaches the soliton solution in the 
continuum limit. As an example, for uj = 0.2, we present in Figure 8, a plot of cos 9i for the 
first three sites (solid curves) and also the corresponding quantities for the sampled soliton 
solution (dashed curves). For uj G (0.424, 1) the (— o+)-brcather continues to the static 
( — h) solution, as demonstrated in Figure 9 for u; = 0.55. These results are summarized in 
the second column of Figure 5. 

It is, of course, no accident that the frequency (a; = 0.424) at which the (o-)-)-breather 
fails to continue to the soliton is precisely that at which the (— o+)-brcathcr begins to 
continue to the soliton. This can be understood by completing the bifurcation pattern of 
the ( — h) solution, using the above results on the continuation of the (— o-|-)-breather, to 
fill in the bottom half of the schematic in Figure 6. As presented in the schematic, there is 
always a bifurcation of the ( — h) solution, but it switches over from the (— o-|-)-breather to 
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the (o+)-breather as the frequency is decreased to the critical value uj = 0.424. At this point 
the (o+)-breather no longer continues to the soliton, but the (— o+)-breather now does. We 
have verified this structure with a number of detailed further calculations, for example, we 
have confirmed the turning back of the continuation of the (— o+)-breather by perturbing 
the (— +) solution as described earlier. 

The lower edge of the band lo G [0.133, 0.424] can also be understood in a similar fashion 

by an analysis of the bifurcation pattern of the ( h) solution, which requires computations 

of the continuation of the ( o+)-breather. These calculations have been performed and 

the results are summarized in the third column of Figure 5. Again there is a soliton band, 
which begins at the frequency u = 0.045, and ends at the start of the (— o+)-breather 

soliton band where to = 0.133. Below the soliton band the ( o+)-breather continues to 

the ( h) solution and above the soliton band it continues to the ( h) solution. 

Computations confirm that the bifurcation pattern of the ( h) solution is analogous 

to that of the ( — h) solution depicted in Figure 6, with an extra — sign inserted into the 
coding sequence of each solution together with a shift in the frequencies. 

Given the above results we are naturally led to the conjecture that the pattern of soliton 
bands continues, with each band covering a smaller range of frequencies, but such that for 
any frequency uj £ (0,1) there is a discrete breather solution whose continuation from a = 
tends to the soliton solution in the continuum limit a ^ oo. 

Generally, as mentioned in the introduction, discrete breathers in networks of coupled 
oscillators fail to continue beyond a certain coupling due to resonance with phonons. In the 
system considered here, the role of phonons is played by spin waves, 

cos{ki — LOpt) 

6ni{t) = I sm{ki — ujpt) | , (4.3) 




which are travelling wave solutions of (2^) linearized about the vacuum n = 63. These 
have dispersion relation 

cjp = 1 + 4asin2(/c/2), (4.4) 

so the spin waves form a frequency band with iVp £ [1,1 + 4a]. Note that the maximal fre- 
quency spin wave {k = tt) is standing, hence of the form (2.4) with frequency ujp = l + 4a. It 



may equally well be regarded as periodic of frequency iVp/n for any n £ Z^. Mathematically, 
a bifurcation in the continuation of breathers may only occur where DF acquires nontriv- 
ial kernel, that is the equation of "motion" ( |2.5D linearized about the breather supports a 
nonzero solution. Of course phonons, being linearized solutions about the vacuum, never lie 
in ker DF strictly speaking. Nevertheless, since the breather approaches the vacuum expo- 
nentially fast as I i| ^ 00, it is generally accepted that the existence of a standing phonon of 
frequency = nuj generically implies nontriviality of ker DF, that is, the phonon is close to 
a tangent vector in ker DF, approaching it asymptotically as \i\ 00. In the present case, 
the standing spin wave lies outside £2, so cannot be close to ker DF, and hence cannot cause 
a bifurcation. In practice this technical point is irrelevant: the numerics are performed on 
a finite chain where such distinctions are impossible to make. Why then do the standing 
spin waves cause no bifurcations in the numerical continuation? The reason is that our 
map F contains explicit parametric dependence on to, so only the standing spin wave of 
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frequency cOp = uj is in the kernel of DF at the vacuum. Higher harmonics, ujp = nuj are 
not linearized solutions of our ansatz derived equation ( |2.5D , and so cannot cause bifurca- 
tions. This is entirely consistent with physical intuition. These breathers, Hke those in the 
DNLS system, are monochromatic: their time dependence contains no higher harmonics so 
one does not expect them to resonate with the spin waves. It is precisely the imposition 
of a monochromatic ansatz ( |2.4D which introduced the parametric uj dependence discussed 
above. 



5 Linear stability 



Linear stabihty of breathers in networks of anharmonic oscillators was studied by Aubry 
in Q. Despite the many superficial differences between spin chains and oscillator networks, 
we will find that the analytic framework developed by Aubry readily adapts to this new 
setting, yielding similar results. In particular, we will prove that (o-l-)-breathers of all 
frequencies oj G (—1,0) U (0,1) must be linearly stable for sufficiently small a provided 
u;~^ ^ Z. Just how small a must be, and the stability properties of the various other 
breather types will be investigated numerically. 

As in [^, it is technically convenient to truncate the lattice to (large but) finite size A'' 
(the results will be independent of A^). Existence of discrete breathers in such a system 
follows from an identical argument to that of theorem |l[ but with £2 replaced by with 
its usual norm. Equation (^|^) along with fixed endpoints no = riiv+i = then define a flow 
on phase space M = (S*^)^. Crucial to the analysis is the fact that this flow is Hamiltonian 
with respect to the natural symplectic structure on M, namely 

n = ^ni (5.1) 

i 

where Qi is the area form on the i-th 2-sphere. The Hamiltonian if : M — > R is the 



truncation of {2A_), again with fixed endpoints. It follows that the flow is symplectomorphic. 

A period T solution of the system is a fixed point of the return map Pt : M ^ M, 
Pt ■ n(0) 1-^ ^{T). Such a solution is said to be linearly stable if the spectrum of its 
associated Floquet map = (d-PT)n(o) • ^n(o)-^ ~^ ^n(T)=n(o)-^ liss within the closed unit 
disk D = {z C : \z\ < 1} . Recall that : 6n{0) 1— > 6n{T) where 6n{t) is the solution of 



the linearization of ( p. 3D about the solution n{t), explicitly, 

6hi = a[ni x (^rij+i + 5ni_i) + 6ni x (rij+i + ni_i)] + {6ni ■ e3)nj x e^. (5.2) 

Since Pt is a symplectomorphism, J^i is a symplectic map, that is 

n{Te6n,J^i6n') = Q{6n,6n'). (5.3) 

It follows that if A S spec J-'i, so are A, 1/A and 1/A. Hence the solution is linearly stable 
if and only if spec J-'i lies on the unit circle dD = {z £ C : \z\ = 1} . 

It is straightforward to construct J^i explicitly for all the breathers considered in sections 
^ and ^ in the uncoupled limit, a = 0. For each nj(0) = (sin^j, 0, cos 9i) define the ordered 
orthonormal basis 

{si = (cos 0i, 0, sin 9i),e2) (5.4) 
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for Tn.(o)'S'^i so that T„(o)Af = 0^ span (e^, 62). Relative to this basis, the symplectic form 
has the usual block matrix expression, namely 

J7 = diag(...,i?(^),i?(|),...) (5.5) 

where R{ip) is the SO{2) matrix which performs a clockwise rotation through angle ■0) 
explicitly, 

RW = ( '^""^ V (5.6) 

y — sm ip cos ip J 

The Floquet matrix of the frequency lo (o+)-breather is 

J'i = diag (. . . , R{T),R{T), S, R{T),R{T), . . .) (5.7) 
where T = 2tt/u} and S is the matrix 

^=((i-VO- 

Similarly, the (— o+)-breather has Floquet matrix 

= diag (. . . , R{T),R{T),S, R{T),S, R{T),R{T), . . .) (5.9) 

and so on. Clearly all these matrices have spec .7-"^ = {e^*^, 1} C dD so the uncoupled 
breathers of all types are linearly stable. 

As a increases from 0, theorem |l] guarantees the existence of a continuous family of 
symplectic maps J^ia (recall that a 1-^ 6*^ is C^) whose eigenvalues depend continuously 
on a. Do these eigenvalues remain on dD? Clearly only coincident pairs may leave dD in 
tandem, by the symmetry properties of specJ-£. Krein theory states that an eigenvalue 
A G dD may bifurcate off the unit circle only if its associated Krein signature is indefinite. 
This signature is defined as follows. Let Ex C C^^ be the A eigenspace of J^£ and Vx := 
Re {Ex e Ei/x) C M^A^. Then on Vx one defines a bilinear form Qx : ® ^ M such that 

Qx{C,v)-=m,:^iv)- (5.10) 

The Krein signature of A is the signature of this bilinear form. Hence, A may bifurcate 
off dD only if Qx is indefinite, that is, the diagonal map ^ 1— > Qx{^,S,) is neither positive 
definite nor negative definite on yx\{0}. 

Applying this theory to (js^) one sees that V^.t = V^-^t = M^^-^, Vi = M^, 

Qe'T(^,e) = -(sinr)c^e, Qi(e,o = (i-^')c?- (5.11) 

Provided sinT ^ 0, i.e. ^ Z, Q(.iT is definite and hence all the eigenvalues e^*"^ must 
remain on dD, at least for sufficiently small a. On the other hand, Qi is only positive 
semi-definite, so no such constraint applies to the double eigenvalue A = 1. Note however 
that for any nonstatic period T solution n{t) of an autonomous dynamical system, n(0) is 
an eigenvector of J-£ with J^£h{0) = n(0) by time translation invariance. It follows that 
one copy of the eigenvalue A = 1 is fixed for all a. The second is also fixed for sufficiently 
small a since to move it would have to (at least) double up due to the symmetry properties 
of specJ^£. Hence, we have proved: 
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Theorem 2 For each N G and for all lu G (—1, 1)\{0} such that lu"^ ^ Z, there exists 
a{u)) > such that for all a G [0,5(0;)) the continued {o+) -breathers, whose existence on 
the N site lattice is guaranteed by theorem are linearly stable. 

Unfortunately, this argument cannot be applied to multisite breathers such as (— o+). 
This is because Vi is now (at least) 4 dimensional and Qi remains only semidefinite. Time 
translation symmetry is not sufficient to fix all 4 copies of A = 1 under perturbation of a, so 

while the a = (— o+)-, ( o+)- etc. breathers are all stable, instability may develop for 

arbitrarily small a. Note, however, that the Krein criterion is necessary but not sufficient: 
it does not guarantee instability of a > multisite breathers, nor of (o+)-breathers with 
uj^^ £ TL. To investigate these issues, we must resort once again to numerical analysis. 

One may vastly simplify the task of numerically constructing the Floquet matrix Ti for 
this problem by transforming to a co-rotating coordinate system. Let 7^ : M ^ 50(3) such 
that 

7^(^) = diag(i2(a;t),l) (5.12) 

where R(t) is the the SO(2) matrix defined in equation (^.6|), and define new variables Uj(t) 
such that 

n,(t) = 7^(^)ui(^). (5.13) 

Then the equations of motion become 

Uj = auj X (ui+i + Ui_i) + (uj • e3)(uj x es) - Awi (5-14) 

where A = {TZ^^'R.){t) = {Tl^^'R.){0) £ so(3) is a constant antisymmetric matrix. The point 
is that all the breather solutions we have been considering are static in this coordinate 
system, namely Uj(t) = (sin6'j, 0, cos ^j), where 6i satisfies ( |2.5| ) as before. Consequently, 
the linearized flow is defined by an autonomous linear system of ODEs, 

5ui = a[ui X (5ui+i + (5uj_i) + 5ui x (uj+i + Uj-i)] + {6ui ■ e3)ui x es - A5ui. (5.15) 

Using the basis {£1,62} for T^^-S^ and writing 

6ui = ai{t)ei + bi{t)e2, (5.16) 

this system reduces to 

:)-(:)Ka- 

where are constant, tridiagonal, symmetric real N x N matrices with components 

Kji^) = a{-5ij+i - 5ij^i + [cos{9i - 9i+i) + cos{9i - ei-i)]5ij} 

+ cos 9i{cos 9i — uj)5ij, 
Kji^) = a{cos(6'i - 9j){5ij+i + 5i,j-i) - [cos(6'i - 6*4+1) + cos(6ii - 9i-i)\5ij] 

+ {uj cos 9i — cos29i)5i J. (5.18) 

The Floquet matrix J^i : 5n(0) 1— > 5n{T) is simply 

J^i = n{T) exp(^A)7^(0)-^ = exp(rA). (5.19) 
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One may easily check that equations (|5.1S| ), ( 5.1S| ) reproduce the a = matrices previously 



obtained. Note that we have not imposed a monochromatic ansatz for the perturbation: 
5ui and hence 5nj are permitted to have arbitrary time dependence. We have merely made 
a particularly convenient choice of coordinates. 

So a period T breather is linearly stable if and only if its corresponding A matrix 
has spec A C iM C C, a criterion which is trivial to test numerically. This should be 
contrasted with the procedure usually employed in numerical linear stability analyses, where 
construction of J-i requires a coupled system of 2N{N + 1) linear and nonlinear ODEs 
to be solved numerically. Of course, the transformation to co-rotating coordinates, and 
the consequent simplification of the Floquet problem, are only possible because we are 
considering anisotropic spin chains which retain a 50(2) isotropy group. The technique 
will not work in chains with fully anisotropic exchange interaction, as are considered in Q. 

The numerical results described below were obtained by computing spec A using various 



standard routines taken from on a 41 site lattice (the results do not differ significantly 
from those obtained with an 81 site lattice). It should be noted that, while the algorithm 
to construct the breather initial profile 6 employed spatial reflexion symmetry, no such 
symmetry is imposed on the spectral problem for A. Hence all possible modes of instability 
are included in the analysis. Given the spectrum {Aj} of A, one constructs 

;u := max(ReAi)^. (5.20) 

i 

The corresponding breather is stable if and only if /i = 0. 

Figure ^ presents a graph of the maximal a{uj) of theorem |2|, that is, the coupling at 
which instability of the frequency lo (o+)-breather first occurs. Three features should be 
noted. The first is that, contrary to expectation, weakly coupled breathers remain stable 
even at resonant frequencies, £ TL. In fact, there seems to be nothing special about 
these frequencies at all from the standpoint either of breather existence (see sections |3| and 
^) or of stability. Second, 5(a;) appears to be finite for all a;: even when lo > 0.424 so that 
the (o+)-breather continues all the way to the continuum soliton (a linearly stable solution 
of the continuum system), it loses stability along the way. Third a{uj) varies very little with 
Lo except for uj close to 1. It appears to grow unbounded as w — > 1, which seems reasonable 
since the uj = 1 (o+)- "breather" is nothing but the trivial vacuum rij = which is stable 
for all a. This limit is numerically inaccessible, however (the breathers tend to spread out 
over the whole lattice as to grows so that finite size effects dominate), so one should treat 
this observation with caution. 



Figure 11 shows plots of fi against coupling a for (o+)-breathers of various frequencies. 
For UJ <C 0.424, the breather remains stable until after it has joined the ( — h) branch, only 
losing stability when this trivial static solution does. After this, /x grows monotonically with 
a. For UJ > 0.424 the breather first loses and then regains stability: a hump of compact 
support appears in the graph of /i. This hump is followed by another of roughly the same 
width but much less tall (note the change in vertical scale), so that stability is once again 
lost and regained. This leads us to conjecture that instability occurs in regularly repeating 
bands, the degree of instability (the height of the humps) decaying exponentially as a 
increases and the breather approaches the continuum soliton. It is interesting to note that 
the transition from stability to instability (or vice versa) generically coincides with a sign 
change in the determinant of DF, where F is the unsymmetrized finite lattice version of 
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Figure 10: The coupling of first instability a as a function of frequency, to for (o+)-breathers. 

the continuation function F defined in (|3.2| ). However no such sign change occurs in the 
Jacobian of the symmetrized continuation function as employed in our continuation 

scheme, so the corresponding eigenvector 66 G ker DF must have odd spatial parity. It 
seems likely, therefore, that stability transitions generically accompany symmetry breaking 
bifurcations in the full (i.e. unsymmetrized) continuation problem. 
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Figure 11: The instability function /u(a), as defined in equation ( [5.20| ) for various (o+)- 
breathers (solid: lo = 0.5; dashed oj = 0.1). The breather is stable if and only if = 0. 
Note that both are stable for small a despite having resonant frequencies. 

Figure ^ shows similar plots for (— o+)-breathers. Recall that no result equivalent to 
theorem § holds in this context, so there is no reason to expect stability even for arbitrarily 
small a. Indeed, we found that for all uj, fi departs from as soon as a does. For lj ^ 
(0.133,0.424) the breather remains unstable until it joins the appropriate trivial branch 

(( h) for u) < 0.133, ( — h) for iv > 0.424), whereupon it becomes briefly stable before 

instability irrevocably sets in. For uj £ (0.133,0.424), the breathers gain stability, then lose 
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it, rather in the fashion of the (o+)-breathers described above. Again, a pattern of repeating 
bands of diminishing instabiUty seems to occur as a grows large and the breather approaches 

the continuum soUton. Plots of /i(a) for ( o+) breathers with uj £ (0.045,0.133) are very 

similar: the breathers start unstable, then gain and lose stability in a repetitive pattern. 
This leads us to conjecture that wherever breathers tend to the continuum soliton as q ^ oo, 
a similar banded stability pattern occurs. 



(a) xio-^ W 




Figure 12: The instability function /x(a) for various (— o+)-breathers (solid: u = 0.3; 
dashed: uj = 0.6; dotted: uj = —0.4). All are unstable for small a, though they have 
nonresonant frequency. 



6 Conclusion 

In this paper we have proved the existence of discrete breathers in ferromagnetic spin 
chains with easy-axis anisotropy, and constructed such breathers by means of numerical 
analysis. The main novelty of our numerical results is that discrete breathers exist inde- 
pendent of the inter-spin coupling a, right up to the continuum limit a oo. This should 
be contrasted with previous studies of spin chains with easy-plane anisotropy jl^, where 
the breathers are again "monochromatic," but the spin waves have a different dispersion 
relation, and apparently do cause trouble at large a. 

The method of proof itself was something of a hybrid: an ansatz was made to reduce 
the equation for breathers from a system of ODEs to a purely algebraic system, but then 
existence of solutions of the latter was proved by continuing a decoupled solution in the 
sequence space £2- Consequently, the result shares nice properties of both the Flach and 
MacKay-Aubry approaches. Since the reduced system depends parametrically on a;, spin 
wave resonances cause no problems, theoretical or numerical (c.f. [^). On the other hand, 
the result generalizes immediately to multi-dimensional spin lattices, as MacKay-Aubry 
style theorems usually do. Whether our numerical results generalize similarly, that is, 
whether a picture similar to Figure H emerges for 2 and 3 dimensional lattices is an interest- 
ing open question. It should be noted that the continuum theory in higher dimensions does 
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support radially symmetric, exponentially localized soliton solutions of the form analogous 
to (|2^) , namely 



n(x,t) = (sin0(|x|) coswt, — sin6'(|x|) sinwt, cos6'(|x|)), (6.1) 



although explicit expressions for 9{r) are not known |11]. 

The ansatz may be interpreted as transforming to co-rotating coordinates, wherein the 
breather solutions are static. This viewpoint greatly simplifies the linear stability analysis 
since one need only solve an autonomous linear system of ODEs. In this way we have made 
an extensive study of the linear stability properties of the breathers. The results suggest 
that those breather families which converge as a — > oo to the continuum soliton experience 
repeating bands of instability of diminishing strength as a increases. We found that the spin 
waves have no influence whatsoever on stability issues, even though the stability analysis 
includes all possible perturbations, including those outside the monochromatic ansatz. Flach 
et al also perform a numerical linear stability analysis for breathers in their spin chains, and 
obtain results broadly similar to ours namely 1-site breathers are found to be stable at 
weak coupling, and certain multisite breathers are found to be unstable at weak coupling. 
The subsequent loss and regain of (in)stability for increasing coupling is not reported - 
perhaps it is special to chains with isotropic exchange interaction. Their survey of the 
breather existence domains is rather less extensive than ours, however, presumably due to 
the much higher computational cost in their systems, so it is possible that similar band 
structures to those found here exist there also. 

One open question immediately arises: that of breather mobility in this system. In 



the continuum limit, there exist travelling solitons which propagate at constant speed |11|. 
Can one find analogous propagating breathers at finite coupling? There seems to be little 
hope of proving anything rigorously about such objects, but one could still hope to study 
breather mobility numerically, focusing on how breather mobility depends on uj and a. 
For a survey of what is known about moving discrete breathers, see the review articles 



[P, 12 1 . One striking observation by Aubry and Cretegny is that breather mobility may be 
associated with certain behaviour of the Floquet matrix Q. In a one-parameter family of 
breathers (such as 6*"), mobility occurs precisely at a value, a' say, where a certain type 
of Krein bifurcation occurs: J^i becomes defective and so-called "marginal modes" appear 
(they span the orthogonal complement of the sum of the J^ia' eigenspaces) . Perturbing 
Oa' in the direction of such a marginal mode produces a slowly translating breather. So 
breather mobility seems to be naturally associated with transitions from linear stability 
to instability. It would be interesting to see whether the stability transitions observed in 
section [s] have such significance. 
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